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We point out that the inclusion of a string component contributing around 5% to the CMB power 
spectrum amplitude on large scales can increase the preferred value of the spectral index ris of 
density fluctuations measured by CMB experiments. While this finding applies to any cosmological 
scenario involving strings, we consider in particular models of supersymmetric hybrid inflation, 
which predict Us 0.98, in tension with the CMB data when strings are not included. Using 
MCMC analysis we constrain the parameter space allowed for F- and D-term inflation. For the 
.F-term model, using minimal supergravity corrections, we find that logK = —2.34 ± 0.38 and M = 
(0.518 ± 0.059) X lO^^GeV. The inclusion of non-minimal supergravity corrections can modify these 
values somewhat. In the corresponding analysis for _D-term inflation, we find log k = — 4.24±0.19 and 
TiFi — (0.245 ±0.031) X 10^''GeV. Under the assumption that these models are correct, these results 
represent precision measurements of important parameters of a Grand Unified Theory. We consider 
the possible uncertainties in our measurements and additional constraints on the scenario from 
the stochastic background of gravitational waves produced by the strings. The best-fitting model 
predicts a _B-mode polarization signal « 0.3/iK rms peaking at ^ « 1000. This is of comparable 
amplitude to the expected signal due to gravitational lensing of the adiabatic i?-mode signal on 
these scales. 

PACS numbers: 11.27-|-d, 98.80-k, 98.80Cq 

I. INTRODUCTION 

The publication of the most recent results from the Wilkinson Microwave Anisotropy Probe Q (WMAP) has 
focused attention in the direction of precision constraints on inflationary models 01 believed to be the origin of the 
initial spectrum of density fluctuations By measuring the power spectrum of the i?-mode polarization of the 
cosmic microwave background (CMB) on large-scales, the WMAP data constrains the optical depth to reionization, 
tr = 0.09 ± 0.03, with consequent improvement to parameters such as the spectral index of density fluctuations, 
Hg = 0.95l!l]Q'Qi[g, which are degenerate with tr. Although ad-hoc inflationary models can be constructed with a 
range of values of Ug, some specific models predict very narrow ranges for making them vulnerable to the very 
tight constraints now available. 

One such class of models is supersymmetric (SUSY) hybrid inflation P|, which has F-term 0,0 and ZJ-tcrm 
variants. These scenarios are particularly attractive since the potential for the inflaton is flat at tree level. It acquires 
corrections from loop effects and when further tree-level, non-renormalisable Planck scale suppressed operators are 
added, and it can naturally meet the slow-roll conditions. Within certain minimal realisations, one is left with 
the choice of only one dimensionless coupling constant k and a mass scale M. One can deduce the amplitude of 
curvature perturbations Pn, their spectral index Us and the dimensionless string tension Gfi from k and M (up to 
a weak dependence on the unknown reheat temperature, Tji). Therefore, these models may be considered as rather 
predictive when compared to other scenarios for inflation. 

One can make an analytic estimate of the scalar spectral index, which yields ris > 1 — l/Nc, where iVo is the number 
of e-foldings (measured from the end of inflation) when a particular observed scale leaves the horizon (or conversely 
when it comes back inside the horizon). For standard estimates of the reheat temperature, iVe « 50, making ris 0.98 
a prediction of the simplest version of these models. If one takes the quoted observational constraint seriously, such a 
model would be excluded at around the 2 — a level under the assumption that only adiabatic density fluctuations are 
created during inflation. However, in these models cosmic strings will be formed at the end of inflation, if the phase 
transition induces the spontaneous breakdown of a U(l) symmetry by the waterfall fields. Since the predicted energy 
scale of these strings is around the grand unified (GUT) scale, they may contribute significantly to the observed 
perturbations Q , creating an interesting phenomenology 0, 0| . 

In this paper, we point out that the inclusion of a sub-dominant string contribution of around 5% to the large 
scale power spectrum amplitude of the CMB can increase the preferred value for the spectral index up to Ug « 0.98 
(and the maximum allowed value at 2 — cr level up to sa 1-02), something which is a generic point valid for all 
models of inflation which produce cosmic strings. Naively, it may seem to be a grotesque violation of Occam's razor 
to have two sources of fluctuations with nearly equal amplitude; in no way do we claim that the data requires the 
additional string contribution in a Bayesian model selection sense. But as we shall describe, a string contribution 
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of the required size arises very naturally in the class of models under consideration here. Note moreover, that these 
models constitute an attempt of fusing together the areas of particle physics and cosmology and may therefore be 
considered to be more attractive when put in the wider context. We then proceed to constrain the parameters k and 
M for specific realizations of i^-term and D-term hybrid inflation models. Our results show that the inclusion of the 
string contribution is critical to determining the correct constraints on the parameters. 

We note that the upper bound on the string tension /i has recently been discussed by a number of authors |llL 
Il4l ITsjl. Their basic conclusion, using a variety of different methods, has been that there is a 2 — cr upper limit 
of G/z < (2 — 3) X 10~^; something which we shall confirm. However, they have ignored the effect strings have on 
the preferred value of Us . Qualitatively similar ideas were pointed out in ref. in the context of constraints on the 
global texture model using the first year WMAP data. At that point in time, the constraint on rig was not as tight 
as is the case now and, therefore, the necessity of including the defect contribution was not so critical. 

II. MODELS AND METHODOLOGY 

In Section III Al and III Bl we introduce the inflationary models we are studying and discuss in particular the various 
contributions to the inflaton potential. In 111 CI we specify the model for the string power spectrum we are using and 
in III DI the details of the Markov-Chain-Monte-Carlo (MCMC) analysis are presented. 



A. F-term inflation 

i^-term inflation |^ is implemented by the superpotential 

W = kS0G-M^), (1) 

where S denotes a gauge-singlet chiral superfield, G belongs to a certain A/'-dimensional representation of the gauge 

group and G is a corresponding conjugate multiplet. In the subsequent discussion, we set A/" = 1, as we are studying 
the gauge group U(l), the breaking of which leads to the production of cosmic strings. 

Having specified the model, we now enumerate the various contributions to the scalar potential, which determine 
the inflationary dynamics. The leading order contribution is the tree-level scalar potential 



\GG - M'^\^ + \SG\^ + \SG\^ 



(2) 



where G, G and S are the scalar components of the respective superfields. While S is to be identified with the infiaton, 
G and G are usually referred to as waterfall fields. 

Infiation takes place along the trajectory where the vacuum expectation value (VEV) of the inflaton obeys S > M 
and where G = G = 0. Vb is completely flat in this direction. However, since Vb = K^Af^, SUSY is broken and 

the mass degeneracy is lifted, such that the superfields G and G encompass mass cigenstatcs of one fermion of mass 
kS, and two pairs of scalars of mass squared k^{S'^ ± A'P) each. This induces the Coleman- Weinberg radiative 
correction |6l[l7| 



(3) 



where Q is a renormalisation scale. For notational convenience, here and in the following it is always understood 
that we take the moduli of complex fields, for example S = \S\, unless explicitly stated otherwise. Added to Vq, the 
Coleman- Weinberg correction lifts the flatness and forces S to slowly roll towards zero. When the critical value S = M 
is reached, the waterfall flelds acquire a negative mass square term which forces them to assume the U(l)-breaking 
VEV G = G = M, while S is driven to zero. For realistic scenarios, the U(l)-group arises at an intermediate stage of 
the breaking of the GUT-symmetry down to the Standard model group Gsm — SU(3)c x SU(2)l x U(l)r- Possible 
candidates are the baryon minus lepton symmetry B — L, the right-handed isospin or the groups U(l)j!f and \J{l)z 
from the embedding SO(IO) D SU(5) x U(l)x 3 SU(3)c x SU(2)l x \J{1)x x U(1)z D Gsm x Z2, see for example 
ref. The spontaneous breakdown of any combination of these U(l)-symmetries at the waterfall transition leads 
to the formation of local (gauged) cosmic strings. 

Another interesting variant is Fjj inflation |l9l l20l | , where the matter of the Minimal Supersymmetric Standard 
Model (MSSM) is not charged under the broken U(l), giving rise to new long lived particles which may loosen the 



gravitino bound discussed below. Moreover, in this scenario the inflaton sector is tied to the MSSM-Higgs and to 
TeV-scale right-handed neutrinos. Note, however, that these additional couplings of the inflaton alter the "minimal" 
form of the Coleman- Weinberg potential JSJ. Nonetheless, the results presented here, in particular that the string 
network allows for a blue spectral index, are also qualitatively applicable to the i^c-model, and a quantitative study 
constraining the model parameters, which might also be accessible at collider experiments, can easily be performed 
using the methods applied here. 

Assuming that supersymmetry is local, supergravity gives rise to corrections of the following form 

^SUGRA = clH^S' + 32n^K^M*-^ + ..., (4) 

where = 87rK^M'*/(3mpj) is the squared Hubble rate during inflation and mpi = 1.22 x lO^^GeV denotes the 
Planck mass. While the term cx 5* is uniquely determined when only allowing for renormalisable terms in the Kahler 
potential, the cj^-tcim can be present for a "non-minimal" Kahler potential. If cj^ is not imposed to be zero, for 
example by some symmetry |2l| , it is expected to be of order one, which may incline the potential to an extent that 
it becomes unsuitable for slow-roll inflation; an observation which is known as the 77-problem [H. t23| . For this reason, 
the minimal SUGRA-case with cj^ ~ is most often considered in the literature. For the studies in this paper, we 
use the potential Q), but note that also the term oc gets a correction factor (1 -I- 0{cjj)), which in turn may again 
be modified by additional nonrenormalisable corrections which are theoretically undetermined. It has been pointed 
out recently that a careful choice of the parameter cjj brings the scalar spectral index into accordance with its 
central value determined from the WMAP3 data for models without cosmic strings. We shall return to this issue 
in section UlIDI 

The curvature-induced correction for SUSY in de Sitter background can be derived to be (2^ 

Since supersymmetry is not protected by non-renormalisation theorems in curved space, the derivative of this con- 
tribution with respect to S depends on the cutoff scale Q^. However, it turns out that for reasonable choices of Q 
between ^/kAI and mpi, the curvature correction may only be significant for large values of k, where the contribution 
of strings to the power spectrum is too large to accord with observation. 

Additional corrections from soft SUSY breaking turn out to be negligible, except for the tadpole term 

Vtp = 2KM^asRe[S] , (6) 

where as should be of TeV-scale. The precise value of this parameter is theoretically undetermined, and moreover, 
its effect also depends on the phase of asS. It turns out that it may become important only for low values, k ^ 
5 X 10-4 [23,123 ■ 

Putting everything together, the full _F-term inflationary potential is given by 

V = Vo + Vcw + VsvGnA + Vr + Vtp ■ (7) 

In the subsequent numerical analysis, we study various scenarios arising as special cases of this generic potential. 

Some analytic understanding of the dependences between the parameters can be gained from the following ap- 
proximation. We introduce the canonically normalised inflaton a by making the choice of phase a = a/2Rc[5'] and 
Im[5'] = 0. A simple estimate for the spectral index Ug can be obtained for large values 10^^) by neglecting all 
corrections other than 

X/ + Vbw = Vb -I- ^n'^M^ In . (8) 



in 
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The number of e-foldings between the time te at the horizon exit of the scale k = 0.05 Mpc ^ and the end of inflation 
tc, where the inflaton reaches the critical value Uc = V^M, triggering the waterfall-phase transition, is given by 



dVy _ 327r-^ 
to '"pi \dcr) ^ K2m2j 



N.= l dtH^-^^M^I dar-) (9) 



where we have used the slow- roll approximation 3Hda/dt = —dV/da and have neglected the lower boundary term 
of the integral. We then find the slow-roll parameter 
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FIG. 1: The computed values of (top-left), logP-R, (top-right) logG/i (bottom-left) and A^e (bottom-right) as a function of 
F-term model parameters k and M for Tr = lO^GeV and g=0.7. 

and consequently, when omitting the contribution from the slow roll parameter e, which is negligible in these models, 



1 - 2r/ = 1 



0.98, 



(11) 



for iVc Si 50. Furthermore, the scale of inflation can now be determined by imposing the observed amplitude of the 
primordial power spectrum, 



,mXi dV/da 



4^ 



^VN^^ w 4.54 X 10"^ . 



V3 



(12) 



from which it follows that the symmetry breaking scale is close to the Grand Unified scale, M « 4 x lO^^GeV. 
For P-ji(k), we have taken here the value determined by the standard six parameter fit (the basic set of four plus 
{log(10i°P-R),nJ), see section HTdI 

We are considering the case where the gauge group is U(l) and therefore cosmic strings are expected to form during 
the waterfall transition. If the string network evolves towards a self-similar scaling regime as observed in simulations, 
then the initial conditions are unimportant and the only important parameter is the dimensionless string tension G/i. 
This can be computed from the model parameters k, and M assuming that the strings do not have superconducting 
currents and results for the Abelian-Higgs model [2^ can be applied. In ref . , it was suggested that one could use 



(13) 



where (3 = /{2g'^), g is the gauge coupling assumed to be 0.7 based on Grand Unification and e{(3) ~ 1.04/3"-^'^^ for 
13 > 10-2 while e(/3) = 2.4/ log(2//3) for f3 < IQ-^. 

In our subsequent analysis, we go beyond these simple estimates, which requires us to compute Ug and Gfi as 
functions of arbitrary k and M. We define P-jz and at fc = 0.05Mpc-^. 
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FIG. 2: The computed values of ria (left) and logPn (right) including the curvature Vr, and tadpole corrections, Vtp, to the 
_F-term model for as = ITeV, Tr = lO^GeV and g = 0.7. 



We parameterize the number of e-foldings at the time when the scale k = 0.05Mpc ^ crosses the horizon, iVe, using 
the unknown reheat temperature which is given by 



2 

3 IQiSGeV 



(14) 



This relation is obtained by assuming that inflation is followed by a matter dominated epoch of coherent oscillations, 
then a radiation dominated epoch with the initial reheat temperature Tr , which lasts until matter-radiation equality, 
and eventually through matter domination until the present epoch. 

The available constraints on Tr arise from the requirements that, within local supersymmetry, gravitinos must not 
be overproduced ^] and that successful baryogenesis takes place. Avoiding the gravitino problem gives 10^° GeV 
as a conservative upper bound for Tr, but depending on the gravitino mass and its branching ratio into hadronic 
decays constraints on the temperature as low as 10^ GeV may be in order. The bound on the reheat temperature 
may be relaxed, when heavy particles decay at late times and release entropy, a feature which is naturally invoked 
within the f n-model Il9l l20| . Scenarios for leptogenesis require lower bounds of > 2 x 10^ GeV for non resonant 
thermal leptogenesis |29l |. 0.3 TeV for resonant thermal leptogenesis |30l | or 3 x 10'' GeV when right handed neutrinos 
are generated directly through the decay of the inflaton [23. Note that in the range 0.3 TeV < Tr < 10^° GeV, the 
number of e-foldings varies just by two. In section [111 Bl we also study the effect a variation of Tr within the above 
bounds has on our results and find it to be small. Throughout the rest of our analysis, we therefore choose the value 
Tr = 10^ GeV, since lower reheat temperatures require unnaturally small couplings of the inflaton-waterfall sector to 
the Minimal Supersymmetric Standard Model (MSSM) matter. 

Including the full potential and performing a numerical study for the minimal SUGRA-case cj^ — and ignoring 
the contributions to the potential Vr and Vtp, it turns out that 0.98 is just a lower bound for the spectral index. 
For large k, the spectrum is turned blue due to the SUGRA-corrections. whereas for small k, the inflaton evolves so 
slow that a quasi scale- invariant spectrum is generated [23, H^, H^, 113, 0, 113 ■ This is illustrated in the Fig. ^ for 
Tr = lO^GeV and g = 0.7. We note that we have used the full formula for Us in terms of the slow-roll parameters rj 
and e, although the contribution from e turns out to be negligible. Moreover, we have computed the running of the 
spectral index, rirun, and the tensor-to-scalar ratio, r; we found that jrii-unl < 10"'^ and r < 10""^ for all the values of k 
and M in Fig.^ We already pointed out that we expect corrections at low k, if we include Vtp, which is illustrated in 
Fig. Elfor as = ITeV. As discussed above, the value of as is not determined theoretically. Qualitatively, one expects 
that for lower as the deviations from the 05 = case only become important for even lower values of k. We have 
checked for example, that the minima of P-jz in Fig. [21 is moved about a factor of two further to the left when choosing 
as = 0.1 TeV. We have also included the curvature correction Vr, which has no significant impact since at large k, 
where it might become important, as the potential is dominated by the minimal SUGRA correction. 



B. D-term inflation 

The D-term model shares many features with T-term hybrid inflation. It is implemented by the superpotential j3| 

W = KsdG, (15) 
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FIG. 3: The computed values of Us (left) and mpi for a fixed Piz, given by 11211 . The solid line has the SUGRA parameter 
g = 10~*, the dotted line g — 10~^ and the short-dashed line g = 10~^. Below g = 10~^ the observable properties are only 
very weakly dependent on g for the range of k considered. The long-dashed line shows g = 0.05, and for small values of k the 
induced values of mpi and g axe increased. 



and by the D-tcrm 



where mpi denotes the Fayet-Iliopoulos mass. These combine to give the tree-level scalar potential 

1 



V 



'GGf + \SG\^ + \SG\'' 



(16) 



(17) 



Note that the ZJ-term is also present in F-term models, but that it just imposes \G\ = \G\ since mpi = 0, while in the 
present case, it is responsible for the spontaneous breakdown of the U(l)-gauge symmetry after inflation which occurs 

at 5^ = 5^ = ^m|j. Above this value for S", G = G = 0, and since Vn = ^ 



pj is non-zero, SUSY and therefore 



the mass-degeneracy between bosons and fermions is broken. The mass eigenstates comprised in G and G are one 

2 

Dirac fermion of mass ks and two pairs of scalars with mass squared k^s^ ± ^mpj, which induce the the radiative 



correction 
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(19) 



Note that a corresponding correction also applies to the Coleman- Weinberg potential in the i^-term model, but it is 
negligible when compared with the minimal SUGRA correction in that case. In turn, the minimal SUGRA correction 
and the possible cj^ term do not occur in the D-tcrm model due to the vanishing of the i<"-terms. When ignoring 
SUGRA or assuming its absence, we can simply set s = S. One can then prove that the VEV of the inflaton at 
horizon exit CTc is proportional to g. For a given P-ji, the induced values of tofi and are then independent of g. 
Numerically, it turns out that the approximation s = S and therefore the degeneracy in the parameter g is good for 
g ^ 0.1 and where the parameter k is within the allowed range, see Fig.|2| 

We note that a tadpole correction does not occur in the D-term model, but there is the curvature term (jSj). Since 
this only contributes significantly for large values of k where there is a large contribution of cosmic strings to the 
CMB-tcmpcrature fluctuations, we do not consider it any further here. 
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FIG. 4: The computed values of Us (top-left), logP-n (top-right), logG^ (bottom- left) and A^e (bottom-right) as a function of 
D-term model parameters k and AIfi for g = 10~^, Tr = lO^GeV and s = S. 

The dimensionless string tension in these models is given by G/i = 27r(TOFi/wpi)^, since the strings satisfy the 
Bogomol'nyi bound. Moreover, the relationship between TVc and Tr is modified to 



iVc « 50 + 



3 IQSGeV 



2 Vff'^Fi 

3 IQiSGeV 



(20) 



As in the F-tcrm case we have computed the important observational quantities ignoring the Vr term, and the 
results are presented in Fig. 0] using g — 10~^, Tr = lO^GeV and s = S. In accordance with the above discussion, 
these results are independent of the parameter g if we ignore SUGRA, and are also a good approximation for the 
minimal SUGRA case when g is small (< 0.1). 



C. Cosmic string power spectrum 



Work in the late 1990s led to strings being excluded as the primary source of cosmic fluctuations |33ll34Ll35ll36| |. since 
they cannot reproduce the observed peak structure of the power spectrum due to decoherence [sj , and the maximum of 
the spectrum is located at the multipole moment £ « 500, since the fluctuations are created at the scale correspon ding 
to the correlation length [s^. More recent work has improved the accuracy of the string spectrum [Tsl Isol l40l| . 
but since the string contribution to the power spectrum is at most 10%, substantially less accurate predictions are 
necessary than for the dominant adiabatic component. 

We have used the model described in ref. ^3 which is an adaptation of that first proposed in ref. [s^- It models 
the string network in the radiation era as a set of line segments with a given length — 0.26 relative to the horizon 
and rms velocity (w^)^ = 0.65, where the functional extrapolation from the radiation era to the matter era can be 
found in ref. ■ The adaptation of ref. also includes the effects of string "wiggles" 0, via the parameter 



/3 = y fj-/T where fl and T are the mass per unit length and tension of the wiggly strings. The value of /3 in the 
radiation and matter eras has been estimated to be /?r ~ 1-9 and /3m w 1.5. The interpolation between the two epochs 
is achieved using the function /3(?7) = 1 + (/?,. — l)a/(77d), where an overdot represents the derivative with respect to 
the conformal time 77. In Fig. |S1 wc plot the angular power spectra for the temperature and polarization predicted 
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by cosmic strings with /3j. = 1.9. This spectrum was computed by averaging over 400 string network reahsations and 
using the WMAP best fit cosmological parameters Q^h'^ = 0.0223, Q^^h'^ = 0.127, h = 0.73 and tr = 0.088. f^b, 
f2c and Jim = + are the densities of baryons, cold dark matter and matter defined relative to critical, and 
h = iyo/(100kmscc^^ Mpc~^). We also show the effect of changing the wiggliness parameter on the temperature 
power spectrum by plotting the ratio of the spectrum compared to (3,- = 1.9 for a variety of values of (3r- We see that 
if 1.3 < pr < 2.8 the effect of changing Pr is at most a 20% effect. We find similar size modifications for sensible 
variations in the other two parameters, and (w^)^^. 
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FIG. 5: Spectrum of anisotropies predicted by the cosmic string model used in this analysis normalized to COBE. The solid 
line is the temperature anisotropy, the dotted line is the iJ-mode polarization and the dashed line if the _B-mode polarization. 
In the right panel we show the ratio of the temperature spectra for j3r = 1 (solid), /3r = 1-3 (dot), (3r = 1.6 (short dash), 
= 2.2 (long dash), /3r ~ 2.5 (dot-short dash), — 2.8 (dot-long dash) compared to /3r = 1-9. 



D. MCMC analysis 

Since computation of the string power spectrum takes ~ 24 hours for each set of cosmological parameters it is 
not feasible to do this at every step of a Markov-Chain-Monte-Carlo (MCMC) analysis. However, since the string 
power spectrum is likely to be only around 5% of the total, if the string spectrum varies by less than, say, 10-20% in 
the range of parameters allowed by the WMAP data assuming only adiabatic perturbations, the error introduced by 
assuming that the string spectrum is unchanged by any variation in the cosmological parameters is less than the 1% 
accuracy claimed by codes such as CMBFAST j4J|| and CAMB This would allow very fast MCMC analysis using a 
single extra parameter Gfi which normalizes the string power spectrum. 

In Fig. we present the ratio of the string spectra for parameters which are 3 — cr away, as defined by the constraints 
on adiabatic models from WMAP, from the fiducial spectrum. We see that there is, at most, a 20% variation in the 
spectrum in each case. Therefore, assuming that the string model itself is correct, it appears safe to ignore the 
variation of cosmological parameters on the string spectrum. 

The MCMC analysis used the May 2006 version of CDSMDMC in order to create chains which were used to estimate 
confidence limits on the cosmological parameters. The basic set of four parameters {ilb^i^, ^^c^i^, tr, ^a}, where 9a 
is defined by the ratio of the sound horizon to the angular diameter distance at the redshift of recombination, were 
used in each case. In addition, we vary the set {log(10^'' P^), tIs, G/i}, or alternatively derive these parameters from 

an infiationary model of interest defined by {M, log k, c|^, mpi, log g, }. In tables ITTl and IIIII of the appendix, there 

is, for comparison, a fit for the parameters {ilhh^ , ^ch^ , tb., a, ^og{10^'^ Pji) , Ug}; we refer to this as the standard six 
parameter fit, there and elsewhere in the text. 

For the analysis of section lTil Al we use three parameters {\og{10^° P-ji) , ng, Gfj.} to describe the power spectrum. For 
the analyses of the F-term and _D-term models, we first tried using the model parameters {M, logK} or {mpijlogK} 
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FIG. 6: The effect of different cosmological parameters on the temperature power spectrum predicted by the fiducial cosmic 
string modeL In each panel we plot the ratio of the spectrum predicted by the the WMAP best fit parameters (with /3r — 1.9) 
compared to upper and lower 3 — cr bounds of these parameters . In the top-left panel, for example, we plot the 3 — cr upper 
bound on ^Ith^ = 0.0244 (solid) and the 3 — cr lower bound on Qf,h^ ~ 0.0196 (dotted). In the top right panel we plot 
n,„/i^ = (0.148,0.100), the bottom-left panel shows h = (0.82,0.61) and the bottom-right panel tr = (0.172,0). 



and computed {\og{10^° P-ji) , ns,Gfi} from them. However, since the range of these parameters allowed by the data is 
a very narrow degenerate hne, the acceptance rate for the Markov Chains was very low. As an alternative we found 
that using log(10^''PT;,) and logK as the parameters allowed for rapid convergence of the chains. The values of M (or 
TOpi), Gfi were then computed as derived parameters. Note that this approach also corresponds to a six parameter 
fit. For the simplest models, one needs one further piece of information, A'o and this was parameterized by the reheat 
temperature, Tr. In addition, where necessary cj^ and log 5 were used as additional parameters. 

We use data from the 3rd year observations from WMAP 0j Oil (WMAP3 and WMAPl are used to refer to 
the 3rd year and 1st year data) and three experiments which observe at substantially higher resolution than possible 
using WMAP. These are the Cosmic Background Imager (CBI) [i^, the ArCminute Bolometer ARray (ACBAR) [s^ 
and BOOMERANG IHIllil. The intrinsic flat priors, listed in Table HJ were chosen to be sufficiently broad to 
incorporate the lines of degeneracy known to exist within the space of parameters. 



III. RESULTS 



In this section, we present the results of the likelihood analyses for various scenarios. The standard six parameter 
fit extended by the string tension Gfj. is discussed in IIII Al the i^-term models with minimal SUGRA in IIII Bl D- 
term inflation in IIII Gl and finally i^-term inflation with non- minimal SUGRA in IIII Dl All results of this section are 
summarized in tables HTl and IIIII in the appendix, so that they can easily be compared. 



10 



Parameter 


Prior 




(0.005, 0.1) 




(0.01, 0.99) 


9a 


(0.5, 10) 




(0.01, 0.9) 




(2.7, 5.0) 


Us 


(0.5, 1.5) 


log K 


(-5.0, -0.3) 


log(TR/10^GeV) 


(-6.0, 1.0) 




(-0.25, 0.03) 


logs' 


(-2.0, 0.0) 



TABLE I: Table of flat priors. The notation (a, b) for a particular parameter gives the lower and upper bounds allowed in the 
fit. 



A. Inclusion of strings allows blue power spectra 
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FIG. 7: Results for a 7 parameter fit to the CMB data. The contours are the 68% and 95% confidence levels in each of the 2D 
planes. 

It was pointed out by the WMAP team that the tight constraint on Ug = 0.951!to oig only apphcs to r = 0, and 
that larger values are allowed for non-zero r. The SUSY hybrid inflation models under consideration here give rise 
to r < 10"'*, but the string contribution to the power spectrum is similar in many ways to that for tensors and our 
investigation of this effect was first motivated by the possibility that there is a degeneracy between Ug and G/i which 
would allow Us ~ 0.98 to be accommodated more comfortably. 

The results of performing a 7-parameter fit (the basic set of four plus {log(10^"P7^), Ug, G/x}) are presented in Fig.|7| 
There is a 2— ct upper bound on Gfi < 3.0x 10~^ (marginalized) which is compatible with previous estimates 
It can be seen clearly that for Gfi < 10^^ lower values of = 0.94 — 0.98 are preferred, whereas for G/i > 10~^ 
values as large as Us — 1.02 are within the 2 — a contour. Larger values of Us also require larger values of Ob/i^; if we 
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fix Ug = 0.98 and perform a 6 parameter fit then G/i = 2.5 x 10~^ is the best fitting value and fib^i^ ~ 0.025. This 
analysis appears to compatible with the intuition described above and suggests that a more detailed analysis of the 
more restricted parameter space for the F- and D-term models is warranted. 

The computed value of — 21og£ = 11302.8 for the likelihood of the 7 parameter fit compares to 11305.5 for a 6 
parameter fit which corresponds to a Ax^ « 2.7. It is clear from this that any sensible Bayesian model selection 
criterion would not favour the 7 parameter model over that with 6 parameters. If we fix ng = 0.98 and fit for Gfi as 
the sixth parameter, then we find that — 21og£ = 11303.0 suggesting that this 6 parameter model gives an equally 
good fit. 

If we impose fib/i^ — 0.020 ± 0.002 as suggested by Big Bang Nucleosynthesis wc find that Gfi < 2.2 x 10"^ 
(marginalized) and Ug = 0.953 ± 0.015. 

B. Minimal F-term models 

The results of an analysis of minimal F-term models in which we have varied k and \og{10^^ P-ji) with = lO^GeV 
and g = 0.7 including the string contribution (with n^, G/i, M and Ng as derived parameters) are presented in Fig.|Sl 
We find that logK = -2.34 ±0.38 and M = (0.518 ± 0.059) x lO^'^Qgy, although there is a strong correlation between 
the two. The best fitting model has — 21og£ = 11303.3, which is close to that of the standard 6 parameter fit. The 
best fit value for the string tension is = 2.6 x 10^^, which corresponds to 6% of the power spectrum amplitude at 
£ = 10. We find that Ug = 0.985 ± 0.004, Gfi = (2.5 ± 0.6) x 10"^ and Q^h^ = 0.0255 ± 0.0009. 




login:] lc)g[)C] log[i<] log[i<] 

FIG. 8: Constraints on the F-term model with minimal SUGRA and g = 0.7, Tr = lO^GeV. On the left are the allowed 
regions when strings are included and on the right when they are not included. In the latter case the value of Gn that would 
be inferred by the preferred values of k and M is presented. It is clear that the inclusion of the string component is critical to 
determining the constraints on k and M. 

We have performed the same analysis but have left out the string contribution to the power spectrum. The results 
of this analysis are also presented in Fig.|Sl The range of values in the k — M plane which are allowed in this (incorrect) 
analysis are much larger and for the best fit, — 21og£ ~ 11308.4 which is significantly worse than that with strings 
included. If strings are included, large values of k are disallowed since they give rise to an excess string component, 
G^i > 3 X lO-'^. 
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We have also performed an analysis with Tr allowed to vary between ITeV and lO^^GeV. The results are not 
changed significantly, which is what should be expected; only modifies iVg logarithmically, and the important 
quantities, Ug and P-ji, are only changed by a few percent. The efi^ect of varying Tr is illustrated in Fig. El 



0,995 




FIG. 9: Constraints on the F-term model with strings (left), without strings (right), with minimal SUGRA and 
allowing Tr to vary between ITeV and 10^"GeV. The coloured dots encode the reheat temperature Tr. 



= 0.7 but 



In addition, we have performed an analysis where the tadpole (with as ~ I TeV) and curvature (with Q — mpi/yoTr) 
corrections to the potential are included. These results are illustrated in Fig.^Jand show a bimodal likelihood surface 
in the k— AI plane. There is an additional allowed region for small k, which is absent in the case without the tadpole 
contribution. This is due to the fact, that for small k, the parameter M and, hence, also the string tension increase 
again, allowing for rig to be close to one. 

Constraints on the model parameters for hybrid infiation have been investigated previously by a number of authors. 
The most straightforward approach is to calculate the value of M for a given k in accordance with the observed 
value for Pn. The string contribution to the temperature fluctuations is proportional to the string tension, which 
is by H13|) a function of k and M. Imposing the reported results (for example, ref. |l3l |') for an upper bound on the 
string contribution therefore allows one to constrain the parameter k. The above approach is valid, since Pn is tightly 
constrained, and for a given k, M may only vary in a small range. Rocher and Sakellariadou js^l obtain k < 9 x 10' 



as an upper bound. However, Jeannerot and Postma |27l,|33 point out that taking account of the corrections (|13|l due 
to deviations from the Bogomol'nyi limit leads to the significant relaxation to k 5; 5 x 10~^, under the assumption 
that strings contribute less than 10% to the power spectrum at ^ = 4; a result which is in accordance with our 
analysis. These approaches do not render a precision determination of the allowed range for the parameter k due to 
the uncertainty in the upper bound on the string contribution, a shortcoming which we resolve here. More important, 
unlike the approach of refs. Is^ Is^ ]. the methods presented here fully address the impact of the presence of the 
string network on the preferred values of the other cosmological parameters, in particular Ug. We note that [2^l33.l&3 | 
appeared before the WMAP3 data, which makes high precision constraints on Ug possible, became available. 

We also note the work by Fraisse 0, i "^1^° has performed a fit for the standard six parameters and the relative 
contribution of topological defects as the seventh, comparable to our analysis in section ITlI AI While the paper based 
on the WMAPl data ^3 contains plots indicating qualitatively the same dependence of il^ and Ug on the defect 
contribution as we present in Fig.[7| such a presentation is not given in the discussion of the analysis of the WMAP3 
data where it is more relevant. Fraisse's study does not, however, exploit the relation between k and the spectral 
index Ug and therefore comes short of giving a precision determination of the parameter k. Moreover, the deviation 
from the Bogomol'nyi limit is not taken into account, leading to a far too tight constraint on k. 



C. D-term models 



We have also investigated constraints on k and mpi for I?-term inflation models, which are quantitatively different 
from those for F-term infiation. We initially fix g = 10"'^ and = lO^GeV. The results including the string 



13 




0.1 ' 1 ' 1 

-5 S -1-5 -3 -1 

loglit] loglur] 




-1 



log[<l log[<l 

FIG. 10: Constraints on the i^-term model with minimal SUGRA, curvature and tadpole corrections with as = ITeV, g = 0.7, 
Tr = lO^GeV and strings. 

component and, for comparison, artificially excluding it are presented in Fig. 1771 

For the scenario with strings, wc find for the parameters logK = —4.24 ± 0.19, mpi — (0.24 ± 0.03) GeV and the 
likelihood — 21og£ — 11305.0 for the best fit model. When compared to the i^-term model, significantly lower values 
for K are preferred, which in turn also induces lower preferred values for mpi in comparison to M. The reason for 
this can be seen when the strings are not included - a wide range of values for k are allowed by the data, but those 
with K > 10~^ correspond to string tensions which are excluded. Since the strings satisfy the Bogomol'nyi limit in 
this case, the string bound is much more restrictive than in the F-term model. When the string contribution to the 
power spectrum is ignored, we find for the likelihood function of the best-fit model — 21og£ = 11307.9. 

Following our discussion in section III Bl the results for g = 10~^ also apply for all values g < 0.1. This degeneracy 
can be seen from Fig|31 For larger values of g, the degeneracy is broken, but these points in parameter space are ruled 
out due to the violation of slow-roll conditions. The 2 — cr upper bound on g < 0.092 when strings are included and 
g < 0.44 when strings are not included; see Fig. 1121 

D. F-term inflation with non-minimal SUGRA 

Here, we study the F-term model with a non-zero parameter cj^ for the SUGRA contribution (0J and we treat 
cjj as an additional free parameter in the MCMC analysis. The results as presented in Figs. 1131 and 1141 show that 
negative values for cj^ lead to an enhancement of the red-tilt of the spectral index, which is why this parameter 
originally was considered When strings are included, it can be seen from Fig. ^| that there is a region for the 
parameter k where a non-zero cj^ is not required. We find the best fit values for the parameters logK = —1.87 ±0.66, 
M = (0.417 ±0.093) X lO^^^GeV, cff = -0.030 ±0.035 and the overall likelihood -2 log/: = 11302.6. For comparison, 
we also show in Fig. 1141 the results for the model without including strings, where indeed negative values for cj^ are 
favoured over the allowed range for n. It is worth pointing out that the allowed range for n reaches up to values 
K w 0.4, whereas in the models with cj^ = 0, wc find k < 0.08. This is because the minimal SUGRA correction, which 
induces a blue tilt of the spectral index, can to some extent be balanced by negative values for cj^. However, for such 
large values of k > 0.1, has to be tuned rather strongly. The poor convergence of the Markov Chains in that region 
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FIG. 11: Constraints on the D-term model with g = 10"'^ and Tr = lO^'GeV. On the left are the allowed regions when strings 
are included and on the right when they are not included. Note the very different scale for G/i with and without strings. 
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FIG. 12: Constraints on the D-term model with Tr = lO^GeV and variable g. On the left are the allowed regions when strings 
are included and on the right when they are not included. 

is apparent by inspecting Figs. [T^ and 1141 

IV. DISCUSSION AND CONCLUSIONS 

We have computed up-to-date precision constraints on the parameters k and Af for _F-term inflation, and k and 
mpi for £)-term inflation. Assuming that the data is correct and that one of these minimal SUSY hybrid inflation 
models is correct, we have shown that one can measure accurately parameters of supersymmetric models at Grand 
Unified scales. In particular, we have shown that the inclusion of the effects of strings is crucial to establish correct 
constraints. For comparison, we have performed most analyses also for the case where the string contribution is 
artificially excluded. Throughout our discussion we have also highlighted potential corrections to the constraints for 
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FIG. 13: Constraints on the F-term model with strings, non-minimal SUGRA, tadpole and curvature corrections, as = ITeV, 
g = 0.7 and Tr = lO^'GeV. 




-3.5 -3 -Z£ -2 -15 - 


-05 









-0.2 -0.1 S -0.1 -0.<35 




-3.5 -3 -2.5 -2 -1.5 


-1 -0.5 








0.92 0.34 0.» C 


as, 1 







FIG. 14: Constraints on the i^-term model without including strings, with non-minimal SUGRA, tadpole and curvature 
corrections, as = ITeV, g = 0.7 and Tr = lO^GeV. 



the most minimal models, in particular we investigated the effect of the leading non-minimal SUGRA correction and 
of the tadpole induced by soft SUSY breaking. 

Without taking account of strings, one may infer from the WMAP3 data Q that the minimal F term models, 
which predict w 0.98, are in tension with the data at 2 — tr level and that the small coupling domain of hybrid 
inflation, where ris w 1 is ruled out at 3 — cr level [3.l55|. The latter constraint would in fact rule out ZJ-term inflation, 
which requires the parameter k to be small in order not to have an excess contribution of strings to the perturbation 
spectrum. However, the self-consistent analysis performed here reveals that D-tcrm inflation is not yet ruled out 
and, moreover, that the minimal six parameter i^-term models with strings fit the data as well as the standard six 
parameter model. 

There are also some theoretical modelling uncertainties associated with establishing the cosmic string spectrum. 
We have shown that if 1.3 < (3r < 2.8 then the effect on the string power spectrum is around 20%, which is ~ 1% on 

the total. Sensible variations of and {v'^)V^ also induce variations of < 20%. It appears that further refinement 
of the string power spectrum beyond this level of understanding is unnecessary for obtaining accurate constraints on 
hybrid inflation models. 

In addition to theoretical uncertainties there are a number of other issues associated with the use of the data which 
need to be carefully assessed: the discrepancy between the data and a SUSY hybrid inflation model with no strings 
and Us ~ 0.98 is only around 2 — a (95% confldcnce level). Chief amongst the uncertainties is how the polarized 
foregrounds are extracted and how line-of-sight effects such as gravitational lensing |5^ and the Sunyaev-Zeldovich 
(SZ) effect are dealt with. For example in the analysis performed by the WMAP team an SZ contribution was included 
in the analysis with an amplitude which was marginalized over. Since this contribution to the power spectrum will 
be increasing with such an analysis will tend to lower the best fit value of Ug. Conversely, taking into account the 
gravitational lensing effect will have a tendency to increase the value of Tig. We have included neither in our analysis, 
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presuming that they will cancel each other out. Once even higher precision data is available, from for example the 
PLANCK satellite, these effects may come to dominate the systematics. 




1 



FIG. 15: Components of the temperature and polarization power spectra for the best-fitting F-term inflation model with an 
artificially large value of r = 0.1 at fc = 0.05Mpc~^. The solid lines are the adiabatic component, the dotted lines the string 
component and the dashed line is the B-mode component due to the gravitational lensing of the _B-mode polarization of the 
adiabatic component. For the adiabatic and string components, the three curves (top to bottom) are the temperature power 
spectrum and that for the E- and B-mode polarization. 

There is an important observational signature of these models which may be in reach in the near future. Inflationary 
models with non-zero r predict that there will be _B-mode polarization on large scales peaking around £ w 100, but 
as wc have already noted the SUSY hybrid inflation models predict r < 10~'^ and a signal this weak is unlikely 
to ever be detected. However, as we have already pointed out the anisotropics created by cosmic strings create B- 
mode polarization since they do not distinguish between scalar, vector and tensor anisotropics. Fig. 1151 shows both 
the adiabatic and cosmic string components to the temperature anisotropics and polarization for a model with an 
artificially large value of r = 0.1 at fc = 0.05Mpc~^. It can be seen that there is a S-mode polarization signal due 
to the cosmic string component, which has an amplitude of « 0.3/xK at around i « 1000. This has very different 
characteristics to that due to adiabatic tensor perturbations and is of similar amplitude to the contribution expected 
due to the conversion of i?-mode polarization into _B-mode by gravitational lensing. In ref. [H^ it was pointed out 
that string tensions of as low as « 10~^ might be detectable in future CMB polarization missions if one is able to 
"clean" the lensing contribution using high resolution observations. 

Finally we should point out that a network of cosmic strings with G/i ~ 10^*" will lead to a number of other 
potentially observable effects. In particular, the decay of cosmic string loops could create a stochastic background of 
gravitational waves (see ref. and references therein) if the dominant decay channel for the strings is gravitation. 
Such a background is constrained by the lack of timing residuals in the observations of pulsars; the most recent 
constraint being that Qgh^ < 2 x 10~^ at frequencies of / ~ 2 x 10~^IIz [H^. If wc assume that the absolute lower 
bound on the string spectrum is given by the "red-noise" spectrum generated by the decay of string loops in the 
radiation era, then one can use various parameters measured by string network simulations [60l . iBlj and the measured 
value of Slni to compute a bound on Gfj, as a function of the loop production size relative to the horizon, a. The 
results of doing this are presented in Fig.^| We see that for small a there is a plateau with Gfj, < 1.4 x 10"'' and for 
larger values of a (which are probably less likely) more tight constraints are possible. These results are considerably 
more stringent than those presented in ref. [s^ since at that time the limit fig/i^ < 9 x 10^^ was used. 

Taken at face value these results appear to further constrain, but do not yet rule out, the i^-term scenarios under 
consideration here to a narrow range of log k w —3 and M w 4 x lO^^GcV, with a qualitatively similar situation in the 
D-tcvm case. There are, however, numerous uncertainties, particularly in the details of string evolution which could 
substantially change the conclusions and therefore at this stage wc feel that it would not be sensible to include these 
observations in our likelihood analysis. Wc have already noted that the anisotropy power spectrum that we would 
observe in the CMB is not that sensitive to the details of string evolution since it is sub-dominant and therefore it 
is unlikely that one would be able to unequivocally rule in or out these models on the basis of CMB measurements. 
Hence, it appears that improved observations of pulsar timing and a significantly better understanding of string 
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network evolution, directed towards the pulsar bound, would be the best way of constraining these scenarios further. 
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FIG. 16: Constraints on G/i from due to the absence of timing noise in pulsars against the loop production size a. The region 
above the solid line is ruled out at 2 — a if Qgh^ < 2 x 10"^ for / « 2 x lO^^Hz. See note added in proof for the explanation 
of the dashed line. 
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Note added in proof 



Ref. |H3 appeared recently claiming that in fact the pulsar constraint is Qgh^ < 2 x 10 ^. If this is the case, the 
constraint on G/i is reduced by a factor of 10 as illustrated by the dashed line in Fig. 1161 



APPENDIX: SUMMARY OF RESULTS 



We summarize here the best-fit parameters together with the 1 — ct confidence intervals for the various i^-term 
scenarios in tabled] and for the D-tcrm model in table IIIll For comparison, we have included in each table the results 
for the standard six parameter fit and for the seven parameter fit including strings. 
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TABLE II: Constraints on _F-term inflationary scenarios. In each case we quote results with the string contribution (s) and without the string contribution to the 
CMB (ns). Model I is the basic fit, with no infiationary input parameters, and is a six parameter (in the case of no strings) and seven parameter (no strings) fit. 
Model II is the basic minimal F-term model, with six input parameters, while model III has the reheat temperature of inflation Tr as an additional parameter. Model 
IV includes curvature and tadpole correction to the potential (and is also a six parameter fit) and model V includes the non-minimal SUGRA parameter cjj as an 
additional parameter. 





Model 


Parameter 


I(s) 


I(ns) 


VI (s) 


Vl(ns) 


VII(s) 


Vll(ns) 




0.0237 ±0.0015 


0.0225 ±0.0007 


0.0261 ± 0.0011 


0.0234 ± 0.0005 


0.0263 ±0.0011 


0.0234 ± 0.0004 




0.103 ±0.007 


0.106 ±0.007 


0.098 ± 0.007 


0.102 ±0.007 


0.098 ± 0.007 


0.102 ±0.007 


9a 


1.044 ± 0.004 


1.042 ±0.003 


1.050 ± 0.003 


1.045 ± 0.003 


1.050 ± 0.003 


1.045 ± 0.003 


TR. 


0.097 ± 0.032 


0.092 ± 0.029 


0.131 ±0.028 


0.119 ±0.029 


0.131 ±0.029 


0.118 ±0.028 




3.00 ± 0.07 


3.03 ± 0.06 


3.04 ± 0.07 


3.08 ± 0.06 


3.03 ± 0.07 


3.08 ± 0.05 


ris 


0.964 ±0.019 


0.956 ±0.016 


0.9993 ± 0.0003 


0.984 ± 0.006 


1.000 ±0.001 


0.985 ± 0.006 


logK 






-4.24 ±0.19 


-2.10 ±0.89 


-4.24 ±0.20 


-1.97 ±0.76 


log(TR/10^GeV) 






0.0 


0.0 


0.0 


0.0 


mpi/lO'^GeV 






0.245 ± 0.031 


0.730 ±0.171 


0.249 ± 0.029 


0.754 ±0.135 


log(3) 






-3.0 


-3.0 


< -1.0 


< -0.36 




< 3.0 




2.56 ±0.61 


23.7 ±8.1 


2.65 ±0.57 


24.7 ± 6.7 


iVo 






47.7 ±0.1 


48.4 ±0.3 


48.9 ±0.3 


49.9 ±0.4 




0.217 ±0.033 


0.237 ±0.031 


0.181 ±0.022 


0.213 ±0.026 


0.179 ±0.022 


0.213 ±0.026 




0.783 ± 0.033 


0.763 ±0.031 


0.819 ±0.022 


0.787 ± 0.026 


0.821 ±0.022 


0.787 ±0.026 


to/Gyr 


13.44 ± 0.26 


13.66 ±0.15 


13.00 ±0.17 


13.45 ±0.10 


12.97 ±0.18 


13.45 ±0.10 




11.2 ± 2.5 


11.3 ± 2.5 


13.0 ± 2.0 


13.1 ± 2.2 


13.0 ± 2.0 


13.1 ± 2.1 


h 


0.77 ± 0.04 


0.74 ± 0.03 


0.83 ± 0.03 


0.77 ±0.03 


0.84 ± 0.03 


0.77 ±0.03 


-2 log £ 


11302.8 


11305.5 


11305.0 


11307.9 


11305.0 


11308.0 



TABLE III: Constraints on D-term inflationary scenarios. In each case we quote results with the string contribution(s) and without the string contribution to the 
CMB (ns). Model I is the basic fit, with no inflationary input parameters, and is a six parameter (in the case of no strings) and seven parameter (no strings) fit. 
Model VI is the basic minimal _D-term model, with six input parameters, while model VII includes the gauge coupling g as an additional parameter, where we give the 
2(j-level upper bounds on g. 
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